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Abstract 

We consider 1/Q corrections to hard processes in QCD where Q is a large mass 
scale, concentrating on shape variables in e + e - annihilation. While the evidence 
for such corrections can be and has been established by means of the renormalon 
technique, theory can be confronted with experiment only after clarifying the 
properties of the corresponding non-perturbative contribution. We list predic¬ 
tions based on the universality of the 1/Q terms, and compare them with the 
existing data. We also identify the scale of the non-perturbative contributions in 
terms of jet masses. 


1 Introduction 


Perturbative QCD constitutes a well defined framework for understanding hard processes, 
i.e. processes characterized by a large mass scale Q (see for example |l|). While the zeroth 
order approximation is provided usually by the parton model and is well defined, the radiative 
corrections can bring logarithmic factors which depend on the infrared cut off and destabilize 
the theoretical predictions. To avoid this problem, one concentrates usually on a set of 
infrared safe quantities which are protected against such contributions. There are powerful 
general theorems on infrared stability of theoretical predictions based on factorization of 
short and large distances, for a review see [j2j. In this note we will consider shape variables in 
e + e“ annihilation which are infrared safe. As an example one may keep in mind the thrust 
T defined as 


T = max , 


E(p»n) 
1 S| Pi | 


( 1 ) 


where p* are the momenta of the particles produced while n is a unit vector. 

While perturbative QCD allows us to calculate thrust as a series in a small expansion 
parameter a s (Q 2 ) 0: 


(1 — T) = 0.335a,(Q 2 ) + 1.02a 2 (Q 2 ) + 


( 2 ) 


analysis of the experimental data at least at moderate Q 2 indicates also the presence of 1 jQ 
corrections !•!• Theoretically, such corrections are exponentially small, 

1 IQ ~ exp(-b 0 /2a s (Q 2 )) , 


where bo is the first coefficient in the /9-function, and are clearly beyond the reach of a purely 
perturbative approach. The only consisitent way to deal with the power corrections f6j] is 
provided by the operator product expansion (OPE). However, this technique applies only to 
a very selective set of physical quantities like the total cross section. 

Nevertheless, most a recently QCD-based phenomenology of the 1 jQ corrections has been 
emerging 0, 1, 0. n 0, 0. It is worth emphasizing that the very existence of the 1/Q 
corrections can be readily understood. Indeed, consider the emission of a soft gluon with 
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energy u ~ A qcd- The corresponding contribution to the thrust is of order: 

Ac 


1 — T, 


r^QCD uj du ,, n 

~ /„ e~“ s(A ' 


QCD 


y QCD 


Q 


( 3 ) 


where the first factor in the integrand comes from the definition of the thrust, du/uj is the 
standard factor for emission of a soft gluon, and the running coupling a s (K 2 Q CD ) is of order 
unit for a soft gluon. Although Tsoft cannot be calculated reliably in perturbation theory, the 
presence of the 1/Q corrections is obvious. In reality, the existence of the 1/Q corrections 
has been established 0, 0, p~4[ . |8| not via such simple estimates but rather by means of 


infrared renormalons |H|. Renormalons, which are a particular set of perturbative graphs, 
allow us to clarify the convergence properties of the expansions like (2) and estimate their 
uncertainty as powers of 1/Q. 

While the conclusion on the very existence of 1/Q corrections seems to be guaranteed, 
to develop a phenomenology of such corrections one needs means to relate them in various 
processes and to fix the overall scale. These are the issues central to the present note, which 


is considered to be a continuation of the analysis started in [11 


It might be worth emphasizing that renormalons per se cannot be a basis for such a 
phenomenology. Indeed, to apply renormalons literally one has to calculate all the terms 
in the expansion (2) until they start to rise in the fashion prescribed by the renormalons. 
The power-like terms appear then as an uncertainty of the asymptotic expansion and are 
dependent on the procedure for subtracting the perturbative contributions. The ambiguities 
of such a procedure are spelled out in Ref |16| and it has never been tried so far. 

The attempts to develop the phenomenology of 1/Q terms are based therefore on a 
mixture of theoretical inputs and heuristic arguments. The first estimation of 1/Q terms 
in shape variables [Q relied on the version of the renormalon technique which replaces the 
renormalons by terms non analytical in the gluon mass squared, p 2 MM- This procedure 
can be thought of only in the lowest order in ot s . This is also true for the modification of the 
technique which assumes freezing of the running coupling at some scale |l(J. The analysis of 
the data indicated that experimentally 1/Q terms are proportional to the p/Q corrections 
found theoretically [§. Another line of development E I is to evaluate the renormalon 
contributions to a given process to all orders and to develop an operator analysis based on 
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an effective theory arising in the eikonal approximation j2|. 

In our previous paper [[□]] we argued that terms of all orders in the large coupling 
a s (h?Q CD ) factorize in the 1/Q corrections to various shape variables into a universal factor. 
We argued furthermore that non-perturbative contributions share this property of factoriza¬ 
tion since they are also associated with distances much larger then 1/Q. As a result, one 
can both substantiate the relations among the 1/Q terms in the shape variables || and try 
to develop a machinery for deriving furthur experimental consequences from QCD. Similar 
relations for soft perturbative parts are derived in Ref [lSj; no experimental consequences 
are claimed however because of the reservations for unknown non-perturbative effects. 

In this paper we first consider the renormalon technique versus the general operator 
product expansion in case of the total cross section when both approaches apply. We check 
on this example the hypotheses concerning the non-perturbative terms made in JTjJ. We give 
a list of experimental consequences which follow from the universality of 1/Q corrections, 
non-perturbative contributions including. Partly the predictions (or their variations) were 
discussed earlier in IBB' Finally, we identify the overall scale of the 1/Q corrections in 
terms of the non-perturbative contribution to jet masses which arise in the Feynman-Field 
type models. The considerations of this paper are restricted to the kinematic configuration 
where the dominant nonperturbative corrections come from the vicinity of the 2-jet limit in 
e + e— annihilation. 

The outline of the paper is as follows. In sect. 2 we review the renormalon technique 
in case of the total cross section. In sect. 3 we turn to an analysis of the 1/Q corrections 
to certain infrared safe observables in e + e _ annihilation into jets and list predictions which 
follow from the universality of the 1/Q corrections. In sect. 5 we compare the predictions 
with existing experimental data. In sect. 6 we establish the correspondence of the renormalon 
based picture for the power corrections with the phenomenological decription of the non- 
perturbative effects in jet physics a la Feynman-Field. 
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2 Renormalons vs. the operator product expansion 


In this section we will try to elucidate which properties of renormalons are shared by the 
non-perturbative conributions. As a test case we choose the total cross section of e + e - 
annihilation since in this case a more general framework based on the OPE is also available. 

The OPE applies in euclidean space-time. The polarization operator in euclidean space, 
P(Q 2 ), is related to an integral over the total cross section cr(s): 

Q 2 f cr(s ) 


P(Q) = 


7T 


Q 7 


( 4 ) 


At large energies the total cross section is calculable perturbatively: 

«s(s) , i 0 on r t a s («) 


cr(s) = (parton model ) 


1 + + i.409(riTZ)2 _ 12 .805(—4—) 3 + .. 


7T 


7T 


7T 


( 5 ) 


where the first three terms in the expansion have been calculated explicitly see, for example, 
the second of reference [jl]. Phenomenologically, P(Q 2 ) can also be analyzed via QCD sum 
rules H which approximate P(Q 2 ) at moderately large Q 2 as 


Q 


2 dP(Q 2 ) 


dQ 2 


(parton model ) • 1 + 


o.(Q 2 ) , ) 2 )n 


7T 


+ 


Q 4 


( 6 ) 


where G“ u is the gluon field strength tensor and the vacuum expectation value of (G 
incorporates non-perturbative effects. The theoretical framework behind (6) is the OPE 
for the T-product of two electromagnetic currents while the main hypothesis is that the 
power like corrections are responsible for violation of the asymptotic freedom once one starts 
to descend from very high Q 2 and to approach Q 2 ~ 1 GeV 2 . The cross section which 
corresponds literally to (S) coincides with ([|) as far as a s correction is concerned and a term 
proportional to 5'(Q 2 ). 

Now, the renormalons fall so to say in between the two representations (5) and (6) where 
the former apparently does not contain power like corrections while the latter emphasizes 
the Q~ 4 terms. Indeed, on one hand, renormalons are particular perturbative graphs with 
n insertions of the vacuum polarization into a gluon line. Therefore, they are included into 
the expansion (5). For large n and a one-term /3-function the renormalon contribution is 
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proportional to 


P(Q 2 ) 


renorm 






g 4 Jo 


crk In 


91 

k 2 




(7) 


where we have retained only the contribution of the extremum of the integrand at 


k 2 eff ~ Q 2 exp(-n/ 2) (8) 

which is reponsible for the n\ growth at large n. 

On the other hand, one can argue that the perturbative calculation is of no relevance 
once n > N cr where the N cr is defined by the condition (fce//)v cr ~ A qcd- Indeed, once we 
approach the Landau pole uonperturbative effects are expected to become important. The 
renormalon contribution at N cr is of order 

P(Q 2 )n„ ~ ^ (9) 

and exhibits in this way the power like corrections, shown in eq. (6). These two facets of 
renormalons reveal mixing of perturbative and non-perturbative contributions so that only 
their sum is uniquely determined, as is emphasized in Refs (R|, [HJ (for discussion see also 

|0i). 

Imagine now that we would like to build up a phenomenology of the power like corrections 
based on renormalons. A straightforward logical way would be to evaluate all the pertur¬ 
bative terms until they start to blow up because of the n\ behaviour (7); define in some 
way the perturbative series and postulate the existence of non-perturbative terms which 
would compensate for the arbitrariness in defining the perturbative part. From the practical 
point of view, such a program would be difficult to implement since, for example, one of he 
most advanced calculations of the perturbative expansion (5) does not reproduce yet the n\ 
behaviour. In principle, such a procedure would reproduce the Q~ 4 behaviour of the non- 
perturbative contribution. The Q~ l corrections would be subordinated however to many 
terms in the perturbative expansion and hence insignificant. The relations among the Q~ 4 
corrections in various channels would be procedure dependent. 

Instead, one can pick up the renormalons from the very beginning by rewriting one-gluon 
exchange with running coupling constant a s [k 2 ) where k is the momentum flowing through 
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the gluon line, as 


1 r 

Q* Jo 


d A ka s (k 2 ) = ^ 


—2cr6o 


da / drk 


A 


QCD , 


(^f)7 


da 


4 — 2a6n 


( 10 ) 


and defining the integral over the pole at a = 2/b 0 as, say, its principal value. Such a 
procedure corresponds to using the expansion (6) along with representing the matrix element 


< a s (G a ) 2 > by its perturbative part |2l 


{a s (G a ^ v ) )pert — 2^2 ^-QCD J 


da 


2t r 2 QCD J 1 — ab 0 /2' 


( 11 ) 


Note that the OPE can be used for the evaluation of the renormalon graphs since the 
momentum flowing through the line is much smaller than Q 2 , see eq. (8). Calculating, 
furthermore, the renormalon contribution to various crosssections one would reproduce the 
relations among the Q~ 4 corrections following from the OPE since the substitution of a 
particular value of the matrix element (11) does not affect the general relations. 

Our next hypothesis would be that non-perturbative contributions are proportional to 
those of renormalons. The reason is that the both are associated with large distances of 
order A qcd and therefore factorize simultaneously. In the case of the total cross section this 
assumption is obviously true since both the soft-perturbative and non-perturbative effects 
are absorbed into the one and the same matrix element < ( G “j,) 2 >. Finally, the whole 
machinery of the sum rules would be reproduced if one assumes that the non-perturbative 
effects enhance the renormalon contributions so that 


<(G“ )V, « <(G“ ) 2 >„ m _ pe „. 


( 12 ) 


The inequality (12) is understood here in the practical sense that one may neglect at mod¬ 
erately large Q 2 high orders in a s {Q 2 ) as compared to the Q~ A terms. This assumption does 
not have a sound theoretical foundation as yet and should be considered as a heuristic one, or 
as motivated by the analysis of the data. It is crucial, however, to build up a phenomenology 
starting from renormalons El- 

Let us also mention two other points concerning renormalons. First, we took it for granted 
that a s in eq. (10) depends on the momentum flowing through the gluon line. This is in 
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fact not obvious and should have been substantiated by a careful analysis. This is a difficult 
point of the renormalon technique. Second, renormalons are pure perturbative constructs 
and they cannot distinguish which hadronic matrix element of (G“j,) 2 is considered. The 
enhancement due to the non-perturbative effects can depend on the hadronic state, on the 
other hand. Consideration of total cross sections is related to the vacuum epectation value 

of (G%,) 2 . 

To summarize, the phenomenological success of the QCD sum rules suggests the va¬ 
lidity of two basic observations concerning non-perturbative terms, namely, that (a) non- 
perturbative terms share the factorization properties exhibited by renormalons and (b) 
non-perturbative contributions are large in the sense that one can match the power-like 
corrections with the first one or two terms of perturbative expansions. 

In the subsequent sections we will try a similar kind of phenomenology on the 1/Q cor¬ 
rections to the shape variables, in which case there is no OPE but the renormalon technique 
is still available. 


3 Universality of the 1/Q corrections to event shapes. 


In this section we summarize and further explore the consequences of the universality of the 
1/Q corrections to event shape variables. 


In |lT]| we argued for the universality of the 1/Q corrections based on the following 
observations: 

(a) To all orders in perturbation theory, the 1/Q corrections in the renormalon technique 
come from soft gluons. 

(b) The soft gluons factorize into a universal factor in the above. Namely, for an observ¬ 
able O we can write: 

(G)i iQ = Ro(a s {Q~))E S oft- (13) 

Thus, Ro are O dependant and perturbatively calculable whereas E so f t are universal and 
are given by 


1 r dk 2 i 


Esoft Q J k\ 


7 eik(ots(k±))k A 


(14) 
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where 7 e ik has been calculated in perturbation theory up to the two loop level [^5], |26|. It 
has also been noted that the same quantity appears in the radiative corrections to many 
hard processes ]27|, [28]. [ 2 ?| 

(c) Combining (a) and (b), together with the discussion of sect. 2, we argue that since 
we are looking at the same e + e _ annihilation process the above is true non-peturbatively as 

well. 

At this point we would like to qualify the above remarks: The considerations of this 
paper as well as the earlier one jlTJ, apply in the kinematic domain near the semi-inclusive 
2 jet limit; i.e. T —> 1. 1/Q corrections, in fact, do not just arise from this domain but 
also from the multi-jet kinematic region. Eq (13) is a statement about such corrections 
as well and emphasizes that multi-jet configurations are associated with powers of a s {Q 2 ). 
The region T —> 1 itself will be sudakov suppressed at high enough energies and for such 
kinematic regimes the multi-jet configurations are important. I 11 the last section of this paper 
we return to a discussion of this point. For the moment, however, we restrict ourselves only 
to the region where the 1/Q corrections near the 2-jet limit are important. 

We shall estimate the numerical value of QE so f t , including non-perturbative terms, in 
the next section and identify the corresponding non-perturbative quantity in terms of jet 
masses in sect. 5. In the remainder of this section we will discuss further physical quan¬ 
tities , observed in e + e~ annihilation which also receive 1/Q corrections and calculate the 
corresponding Ro- 

We turn first to a discussion of the jet opening angle. The energy weighted jet opening 
angle is an infrared safe quantity and is defined as 


(15) 


/ . 2 \ / VEisin Si . 

{sin 77 ) = (---) 


where E t and Si are the energy and the angle with respect to the jet axis of i-th particle, 
respectively. Here we would like to demonstrate the existence of 1/Q corrections to this 
observable using the renormalon technique. 

Consider the contribution to this quantity from the gluon emission process the crossection 
for which is:( k is the gluon momentum and pi and P 2 those of the quark and antiquark 














respectively) 


da 


C f 


~CXa 


/y>2 I ry> 2 
JU->Xj 2 


(16) 


a 0 dxidx 2 2n (1 — oq)(l — x 2 ) 
where ay = 2 Ei/Q. Let us change variables from x\,x 2 to x 2 and x = k^/Q More explicitly: 


x = 


y^l - xO(l - x 2 )(xi + 0 ) 2 - 1 ) 


3)2 


(17) 


supposing that the gluon is emitted into the same hemisphere as the antiquark with energy 
E 2 . Then: 

1 rtn ( t 'r ~ \ (nr-. -I- — 1 'l 1 

(18) 


(xl +x%)(xi + x 2 - 1) 1 

"Of c 


with 


and, 


no dx 2 dx it (1 — x 2 )(x 2 — 2(1 — aq)) x 

X\ — 1 


X2 X2 
2 + 2 ' 


4a: 2 


1-0)2 


0 < x < 


v/l - r- A 


(19) 


( 20 ) 


Anticipating that the 1/Q correction comes from the region of soft gluons ((x 2 , X\ ~ 1, x% —> 
0) we see that for gluons Esin 2 5 is 0(to) whereas for the energetic quantities it is O(co 2 /Q). 
Hence we keep only the contribution of gluons as far as the 1/Q terms are concerned. 

Thus, 


i Esin 2 5 

K Q 


) = 


dxo 


rl/ 271 ^ dx 2 C F 


-ot s 


x 2 (xf + xf)(xi + x 2 — 1) 


( 21 ) 


Now to get the renormalon contribution, we led a 
x 1 , x 2 ~ 1 and x small. Then: 


x n (1 — x 2 )(x 2 — 2(1 — xi))(2 - aq - x 2 ) 

> a s (kj_) [BT] and consider the limit 


, Esin 2 5 

1 Q 


AC* 


dx 2 




dx ■ x- 


Oi s (k 2 


7T JO 




4x 2 


-X2 


( 22 ) 


Interchanging the order of integration we get: 


, Esin 2 5 

‘ Q 


) « / dx ■ x 

Jo Jo 


1—Ax 2 


dx- 


4 


a s {k j 


7 r (1 — x 2 ) 2 + x 2 


7r rQ/ 2 

—— / dk\ 
2Q Jo 


4CV 


7T 


« s (fc ± ). (23) 


The general features at the origin of these 1/Q corrections are the same as discussed in the 


previous paper fly]], i.e. they come from a region of soft gluons and close to the 2-jet limit. 
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As final example of the 1/Q corrections to jet properties let us consider the expectation 
value of the fractional energy 77 = 2 E(5)/Q emitted inside a double cone of half opening 


angle S centered around the thrust axis. This too is an infrared safe quantity |24|. We 
consider the n-th moment of 77 : 


A/"(Vi) = J 


(24) 


where p( 77 ) is the probability density describing the emission of the fractional energy 77 inside 
the cone. 

We write this moment as 


A/”[7) j = 1 - /(I - (r/ tdri 


(25) 


where the integration region is one with at least one parton outside the cone. Thus for a qqg 
final state: 

^~i 2 2 

M n (5) = 1 - -^-a s /(l- 7f) — dxidx 2 . (26) 

We are interested in determining if there are any power corrections linear in Q to this 
quantity. For this, we again concentrate on the region of soft gluons. We note that in the 
above the integration region is subject to the following constraints for the gluon angle and 
energy: 

77 = 1 if sinO < sinS; r) = 1 — x 3 if sinO > sin5. (27) 

For soft gluons we expand (1 — rj n ) ~ nx 3 and thus we are let to consider: 

ry> 2 I rp 2 

Jb “I Jb 2 


nC F - - 2 1 - 2 


2 tt 


-a s J dx 2 dx 1 —— 


(1 - aq)(l - x 2 ) 




subject to the constraint (]27|) . In the soft gluon region this becomes in terms of the variables 
x 2 ,x: 

2Cf 


(M n (5)) i/Q = -n—^J dx 2 J 

The integration region is determined from: 


dx a s (k |) r x 2 , 

sV ±J {l-x 2 + - -— 


x 1 — x 2 


1 - av 


^(1 -x 2 )sinS < x < i y/l - x 2 ; 0 < x 2 < 1 . 


(28) 


(29) 
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We next interchange the order of integration and concentrate on the small x region to extract 
consistently the 1/Q term using the renormalon technique. We get: 


r-ksinS rjrp fl—Ax 2 

(M n (S)) 1/Q « — n ~T~ ( - / „ <*s(k] 


n Jo 


X J 1—1 — X 2 

sind 


1 x 

{1-X 2 + 


4 C F 


•Q-sinS 


= — n- 


dk±a s (k ]_) — n 


2C F sind r% sinS 


7 r Qsind Jo 7 T Q Jo 

The conditions for the applicability of this calculation is: 

hAqcd 


1 — x 2 

dk±a s (k 2 L ) 


Q ■ sind 


< 1 . 


(30) 


(31) 


Thus we expect non-perturbative corrections to the quantity M n (S) which for small 6 goes 
like n/Q5. Corrections of the form 1/Q5 for small 5 in jet processes was discusssed in ||. 

We should again emphasize here that in order to compare these predictions with the 
data one has to make sure that the kinematic region under study is such that the 2-jet 
configurations are not strongly sudakov suppressed. 


4 Comparison with experimental data. 


In this section we compare the theoretical predictions based on the universality of the 1/Q 
terms with experimental data. 

At this point, we have four observables (1 — T), ( C), (aQ and (Esin 2 5/Q) on which the 
idea of the universality of the 1/Q corrections can be tested experimentally. The results are 
shown in Table 1 below. In listing the predictions we have used the following observation 
consistent with our approach to renormalon phenomenology discussed in section 2. In cal¬ 
culating < 1 — T >, we proceed as in |[I| and obtain the corresponding 1/Q term in the 
lowest non trivial order. To this order this is the same as the 1/Q correction to the heavy jet 
mass, M/. Since the effect comes from the soft region where the coupling is order unit, we 
must resum all such contributions. In higher orders, the light jet mass Mf will similarly give 
a non vanishing 1/Q correction. Thus we must include this in our non perturbative estimate 
of < 1 — T > as well. Consistent with our discussion in section 2 we must therefore allow for 
a non perturbative enhancement factor of 2. This is taken into account in our predictions 
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for the non perturbative contributions to < 1 — T > as well as to < C >, and to< ol/vt >• 
We shall come back to a discussion of this point in sections 4 and 5. 

Table 1 


Observable 

a s expansion 

1/Q correction {GeV) 

R 0 

QE soft {GeV) 

(1 -T) 

0.335a + 1.02a 2 

~ 1 

2 

~ 0.5 

(C) 

1.374a s + 4a 2 

~ 5 

37T 

~ 0.53 

<J r p 

a s / 7 r 

~ 0.8 

71 

2 

~ 0.51 

Esin 2 5 N 

( Q } 

2 a s / 7 r 

1.20 ±0.25 

71 

0.38 ±0.08 


The estimate of the 1/Q correction for the first three observables is taken from ||. This paper 
contains also comparison of theoretical predictions, in a different form, with experimental 
data treating one of variables as an input. In the next section we shall be able to fix the 
overall scale as well. Data on the energy weighted jet angle is from [[|. The energy range for 
the last estimate is also different from the first three: 8-30 GeV for the latter compared to 
upto LEP energies for the former. We should note that using the results of |3j] for the 1/Q 
correction to thrust, we would find for the quantity QE so f t = 0.30 ± 0.08. 

We did not list prediction ( PH ) because it cannot be confronted directly with experimental 
fits to 1/Q corrections. However, the data on the moments exist for various values of S, n [Q . 
An inspection of the corresponding experimental curves reveals that the data do not show 
large non perturbtaive corrections for small 5 as predicted in (|31|). Taken at face value, the 
data are the most serious disagreement with theory that we End. Updating of the analysis 
of the data appears desirable. 

5 The overall scale of the 1/Q corrections. 

In this section we identify the non perturbative scale associated with the 1/Q corrections 
in e + e _ annihilation, i.e we infer the non perturbative counterpart of QE so f t . In order 
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to motivate this identification, we briefly recall the discussion of thrust from |l^]. There 
it was shown that to all orders the 1/Q corrections come from soft gluons and from the 
neighbourhood of the two jet limit. In particular, we identified from simple kinematics and 
from the infrared safety of T, that to all orders, as far as the 1/Q corrections are concerned: 


Ml 


<l-r>VQ = f + 


(32) 


Here we consider two hemispheres divided by a plane perpendicular to the thrust axis and 
M'l and M| are the invariant masses in these. In the limit of soft gluons Mi = Y/j^Pikj, 
with k the momenta of the gluons.Thus in perturbation theory, (1 —T)\/q is essentially twice 
the jet mass squared. By using the factorization of soft gluons, we have upto corrections of 
0(a s (Q 2 )): 

{l-T) 1/Q = 2E soft (33) 

with, 


QE sof t | 


renormalon 


dkj 


■'feik{oi(kl))k A 


(34) 


where, in the integrand we pick up the contribution of the renormalon pole in a s . 

It might worth emphasizing that, once it is established that the 1/Q corrections are 
associated with soft gluons the relation (|33|) between the jet mass and < 1 — T > becomes 
purely kinematical. It is inherent to hadronization models (for review see Bi). or to soft 
perturbative contributions (see, e.g., |29|| ). 

Now, we use this kinematical nature of relation (1331) to fix the overall scale of non- 


perturbative contribution to the 1/Q terms. Namely, we argued in [|lT|] and in the earlier 
sections that the relation which are true to all orders in a s normalized at low mass scale of 
~ A qcd should be true non perturbatively as well. Thus, we are naturally led to the non 
perturbative identification of Q 2 E so f t , near the 2-jet limit, as the average non perturbative 
correction to the square of the jet masses. Phenomenologically the hadronization correction 
to the jet masses is parametrized thus: 


(MU = A Q 


(35) 
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which corresponds to jet momenta receiving a negative correction of order A. For example, 
in the ’’tube” model of hadronization one finds ( see 0 for a review), 

A = j <Pp±p(p±)p± (36) 

where p(p±) gives the pj_ distribution of hadrons in a rapidity-pj_ ’’tube”. In particular, 
according to the analysis quoted in ||, from the experimental data on jet masses, one finds: 

A ~ 0.5 Gev. (37) 


In fact in simple physical models one can directly see the non perturbative connection be¬ 
tween the l/Q corrections to the observables and to the quantity A. Consider for instance the 
energy weighted jet opening angle ( sin 2 r /). Let p(z,p±) denote the appropriately normalized 
distribution of hadrons in a jet with longitudinal momentum fraction z and perpendicular 
component p±. Thus, if ps is the component of the hadron momenta along the jet axis, then 
z = 2 pzIQ. Whence, 

(,my> = f/dzj rfWC.Px) + pi (38) 

The l/Q terms can only come from the neighbourhood of z = 0, hence: 

(sin 2 p) = / dW(p.i) [ d \/J^ +A < 39 ) 


where, p(0,pj_) = p(p±). Thus, 

(sin 2 p) 



d 2 p±p{p±)p± 



(40) 


The closeness of the two approaches is hard to miss. The reason for this similarity is that 
both renormalons and the standard hadronization picture identify bounded intrinsic p± as 
a manifestation of the non perturbative effects. The corresponding derivations of the l/Q 
corrections therefore become identical in this approximation, with the replacement ^ e ik/k\ —> 
7 rp.( Compare equations (22) and (39)). 

As a final example consider the moments M n {5) discussed in section 3. The (sM)' 1 
dependance is easily understood in the hadronization picture. Indeed, hadrons with energies 
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upto p±/sinS leave the cone with opening angle S and reduce the flow of energy inside this 
cone, in agreement with the renormalon picture given in section 3. ( See equation (30)). 

To summarize, we have proposed to identify the non perturbative scale from pure kinemti- 
cal considerations and the theoretical input on the dominant contribution of soft gluons to 
1/Q corrections. So far the absolute scale of non-perturbative terms was fixed only within 
the operator product expansion in the euclidean region |J.The generalization of this idea 
to Minkowskian space has thus far not produced a handle on the phenomenology of 1/Q 
corrections since the operators become nonlocal || [E|. It is amusing that just the same 
transition from euclidean to minkowskian space which is responsible for the change of local 
operators into nonlocal ones brings new tools based on the minkowskian kinematics. How¬ 
ever, as emphasized earlier,our identification can be correct only to order 0(a s (Q 2 )). In the 
absence of the OPE the generalization to higher orders in a s (Q 2 ) is not straightforward at 
all, as we discuss in the next section. 


6 Comparison with other methods. Conclusions. 

In this paper we have attempted to bridge the renormalon picture of 1/Q corrections with 
the phenomenological hadronization models. The considerations of this and the previous 
paper |TT| may be thought of as providing a field theoretic argument for the justification, 
near the 2-jet limit, of these hadronization models which were so far considered at a purely 
phenomenological level. Conversely, one may turn the argument around and use the success 
of the phenomenological models to justify the picture of renormalon inspired phenomenology 
discussed in section 2. 

It might be worth mentioning that not all renormalon-based approaches end up with the 
same kind of phenomenology as outlined above Q. Indeed, consider models with (fictitious) 
gluon mass // || or with freezing of the coupling constant at a scale pi H- The results based 


on this model were reviewed recently in |3(J. It is no surprise that many predictions coincide 
with what we have. A more detailed consideration unveils however that this coincidence is 
rather formal. The easiest way to appreciate the difference is to consider predictions for the 


1 The following discussion is added following the suggestion of a referee report. 
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difference of masses of a heavy and light jets. According to |8|, |10|, |30 

/i rp\ _ 16 a s /i 

~ \ L 1 h iq ~ 

i IQ 


Ml 


Ml 


) 


= <^> 

K Q 2 ' 


3tt Q = ^ ^ 


' Q 2 HQ 

Note that the prediction for (1 — T)i/q coincides with ([33]) upon identification of our QE so f t 
with (8 /3n)/iia 0 (ni). This identification works in some other cases as well. However, the 
central point point about eq (|4T| ) is that adding or subtracting Mf does not change anything. 
This(M 2 = 0) is inherent to models with a gluon mass or with the frozen coupling [[3D] . 
The reason is that these models generalize to non-perturbative physics basing on one-loop 
calculations. 

On the other hand, within the approach developed in this paper we reproduce in the 
leading approximation the standard hadronozation picture as described, say, by the tube 
model. Therefore, the {Ml — Mf) difference depends on the distribution of nonperturbative 
momenta and vanishes if this distribution is dominated by a characteristic value. It has been 
known for some time 0 that this picture is consistent experimentally: 


Ml-M 2 


) < 

HQ 


' Q 2 


Q 2 1 IQ V l/Q 

Similarly, our prediction for the l/Q corrections reads as 


(42) 


U-^l/Q 


Ml + Ml , 
' Q 2 1 


i IQ 


2-M 

v ; 


i/Q 


(43) 


This difference of a factor of 2 compared to ([IT]) is not immediately manifest because of 
arbitrariness in normalization of fij. 

One might wonder how it is possible to have different pictures within the one and the same 
renormalon approach. The point is that, as emphasized in sect. 2, the enhancement factor 
for nonperturbative effects defies naive analysis of divergencies of perturbative expansions. 
In particular, the model with gluon mass makes non perturbative effects, as imitated by 
non analytical terms in mass squared, subordinate to the perturbative contributions in the 
most transparent way. Indeed, since the ji/Q correction is associated with a small fraction 
of the phase space this term is always a minor effect on the background of the ordinary 
perturbative contribution. Within this model, it is inconsistent to consider both Mh and Mi 
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as nonvanishing unless at least a 2 corrections are included. In fact, more careful analysis 
would reveal that even higher orders are needed for consistent treatment of 1/Q terms. 

Within the aproach which allows for enhancement of non perturbative effects there is 
no difficulty to retain 1/Q terms even when the perturbative terms are small numerically. 
Similarly, in QCD sum rules the condensate terms dominate over pure perturbative contri¬ 
butions at moderate Q 2 (see sect. 2). It might be worth noting that according to Ref the 
corrections to the prediction of the model with the frozen coupling can be as large as 50%. In 
the renormalon language, our eq (|33|) corresponds to summation of all these large corrections 
so that the remaining uncertainty is of order a s (Q 2 ) [pd], [L2|]. In a sense this corresponds 
to summing the soft contributions from the multi-jet configurations as well. However, this 
summation cannot be done in the model with frozen coupling since the freezing has no clear 
field theoretical realization. 

Let us also mention paper [[nj which tests, in particular,the conclusions of references 
11, [| on the example of large nj. Although an agreement is found in lowest orders eval¬ 
uated explicitly, the question is raised on possible existence of a s (Q 2 ) • InQ 2 corrections in 
higher orders, due to multijet events. It might be worth emphasizing, therefore, that all the 
statements on two-jet dominance which are made above are to be qualified for the energy 
being not so high as to make the effect of the Sudakov suppression appreciable. If not,then 
terms like InQ 2 ■ a s (Q 2 ) do appear. Moreover, these terms cannot be dealt directly within 
the renormalon technique. The way out is well known however. Namely, one considers the 
resummed cross section and renormalon corrections to it 0 so, mi- For example, in case 
of thrust instead of < 1 — T >i/q one has to consider 


(e 


— i/(l—T) 


i IQ 


(44) 


where the parameter v is large enough to ensure the dominance of the two-jet configuration 
despite the Sudakov suppression. In our discussion of the experimental data we did not 
introduce v since the fits to 1/Q corrections [§, |10j] do not indicate any energy dependence. 

Another issue raised in Ref ]3T[ is the search for a proper choice of a shape variable which 
would be free of 1/Q corrections and hence could be used for reliable extractions of the value 
of a s . We would like to notice that consideration of spherocity could be promising from this 
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point of view. The spherocity S is defined as |32| 


S = (- 

,7T 


mm 


^i\PT 

Si I Pi I 


Perturbatively, it is given as [33 


_ 64 a s 

/ 9 o 

7T Z 07T 


229 3 

-h 64/n- 

9 2 


(45) 


(46) 


Explicit calculation shows that the spherocity is free from 1/Q correction in the leading 
order in a s . Thus to this accuracy: 


CXc 


(S) = Ji(3.28) + 0{1/Q 2 ) 
3ir 


(47) 


fnstread of giving details of derivation let us remind the reader of the simple estimate of the 
non perturbative contribution to the spherocity [[33| : 

'4\ 2 < p T > 2 


A 

non— 


non—pert 




< n(Q ) >' 


InQ 2 


(48) 


Q 2 ; Q 2 

where n(Q) is the multiplicity and < px > is the characteristic non perturbative transverse 
momentum. Since we know that renormalons do reproduce the basic features of the phe¬ 
nomenological hadronization models, it is no surprise that that there is no 1/Q renormalon 
here. Moreover, this kind of argument appears to work for multijet events. Thus, there are 
good chances that the spherocity is indeed free from 1/Q corrections to higher orders in 
a s (Q) as well. 

As we emphasized above, all of our considerations until now have only dealt with the 
leading 1/Q corrections to the quark jets. For gluon jets, the corresponding anomalous 
dimension is different in perturbation theory. Thus it appears likely that a different 
non perturbative parameter will be needed to describe situations dealing with gluon jets. 
Continuing in a similar spirit,we recall that our estimates of the 1/Q corrections have been 
carried out in the kinematic region where the dominant contribution comes from the neigh¬ 
bourhood of the 2-jet limit. However it is clear that though suppressed, three jet processes 
will have to be taken into account in the next order of approximation,i.e order a s (Q 2 ). New 
non perturbative parameters could arise here. It is a challenge to implement such corrections 
into the developed framework. 
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